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Abstract 

In order to estimate the effects of local structure on the Hubble parameter we calculate the low- 
redshift expansion for H(z) and ^ for an observer at the center of a spherically symmetric matter 
distribution in presence of a cosmological constant. We then test the accuracy of the formulae com¬ 
paring them with fully relativistic non perturbative numerical calculations for different cases for the 
density profile. The low red-shift expansion we obtain gives results more precise than perturbation 
theory since is based on the use of an exact solution of Einstein’s field equations. For larger density 
contrasts the low red-shift formulae accuracy improves respect to the perturbation theory accuracy 
because the latter is based on the assumption of a small density contrast, while the former does not 
rely on such assumption. 

The formulae can be used to take into account the effects on the Hubble expansion parameter due 
to the monopole component of the local structure. If the H (2) observations will show deviations from 
the A CDM prediction compatible with the formulae we have derived, this could be considered an 
independent evidence of the existence of a local inhomogeneity, and the formulae could be used to 
determine the characteristics of this local structure. 
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I. INTRODUCTION 


The standard cosmological model is based on the assumption that the Universe is homoge¬ 
neous and isotropic on sufficiently large scales. Nevertheless local observations could be strongly 
affected by local structure as shown for example in [l], and it is important to study its effects. 
The analysis of luminosity density data [ 2 I has provided some strong experimental evidence sup¬ 
porting the existence of local inhomogeneities, but it would be important to confirm it using 
another observable such as the baryonic acoustic oscillations (BAO) measurements 3-8]. The 
BAO scale allows in fact to determine the expansion rate of the Universe H(z) independently 
from the luminosity distance and as such provides an important source of information about our 
Universe. If the H(z) estimations obtained from BAO observations data will show deviations 
from the A CDM predictions this could be considered an independent evidence of the existence 
of local inhomogeneities. This motivates the calculation of a low red-shift formula for H(z), 
able to take into account the effects of inhomogeneities which cannot be fully modeled with 
perturbation theory, as some of the inhomogeneities found for example in sj. A low red-shift 
expansion based on the use of exact solutions of Einstein’s equations is in fact valid also for 
large values of the density contrast or of the gravitational potential. 

The effects of a local inhomogeneity on cosmological observations have been studied already 
Jdll such as equation of state of dark energy or the luminosity distance 
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24], It has been shown for example that the value of the cosmological constant could 
be affected significantly by the presence of local inhomogeneity seeded by primordial curvature 
perturbations fj, which could also lead to the wrong conclusion of a varying equation of state 
for dark energy while only a cosmological constant is present 19|. The origin of these effects is 
that spatial inhomogeneities can change the energy of the propagating photons, contaminating 
the cosmological red-shift due to the universe expansion, and consequently introducing some 
errors in the estimation of parameters based on cosmological models which ignore the effects 
of the inhomogeneities. 

In this paper we will focus on the low red-shift effects of inhomogeneities on the Hubble 
expansion parameter, adopting an analytical approach based on the use of an exact solution of 
Einstein’s equations to model the local structure. We first derive the red-shift expansion of the 
geodesics equations, and use it to obtain the expansion of H(z). We then compute a formula 
for the relative difference between the A CDM and the inhomogeneous case. Finally we 
compare the formulae with the numerical calculations based on the integration of the Einstein’s 
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equation and the geodesics equations, finding a good agreement. We also check that the low 
red-shift expansion formulae are more precise than the perturbative calculation, especially when 
the density contrast is larger. 


II. MODELING THE LOCAL UNIVERSE 


We use the LTB solution to model the monopolc component of the local structure 32 

,2 , (R'(t,r)) 2 dr 2 


ds = —dt + 


+ R(t, rydfl 2 , 


( 1 ) 


1 + 2 E(r) 

where R is a function of the time coordinate t and the radial coordinate r, E(r ) is an arbitrary 
function of r, and R'(t,r ) = d r R(t,r). The Einstein’s equations imply that 
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2 E(r) | 2 M(r) | A 
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P(t,r ) = 


2 M' 
R 2 R' ’ 


( 3 ) 


where M(r) is an arbitrary function of r, R = dtR(t, r ) and we choose a system of units in 
which c = 8nG = 1. 

To compute H(z) we need to solve the radial null geodesics j^l 


dr yj 1 + 2 E(r(z)) 
dz (1 + z)R'[t(z),r(z)} ’ 
dt R'[t(z),r(z)] 

{1 +z)R , \t(z),r{z)\'‘ 


and then we substitute in the formula for the Hubble parameter in a LTB space 38, (39] 


( 4 ) 

( 5 ) 


where 


H(t,r) = ^Hj_(t,r) + ^H\\(t,r), 

( 6 ) 

H(z) = H(t(z),r(z )), 

( 7 ) 

TT (f r] - 

( 8 ) 

H,tt r)- k ' (t ’ r) 

111 , j “ R'(t, r) ‘ 

( 9 ) 
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The analytical solution can be derived 40, [41] introducing a new coordinate rj = r)(t,r ), and 
new functions po( r ) and k(r) given by 


dr\ 

dt' r = 
Po(r) = 

k(r ) = 


r 1 

R = a ’ 
6M (r) 

^3 ’ 

2E(r) 


( 10 ) 

( 11 ) 

( 12 ) 


We will adopt, without loss of generality, the coordinate system in which p 0 (r ) is a constant, 
the so called FLRW gauge. We can then express Eq. (2) in the form 


Gy) =-*(->’ +fa+ £a‘, 


(13) 


The coordinate r], which can be considered a generalization of the conformal time in a homo¬ 
geneous FLRW universe, is defined implicitly by Eq. (10). The relation between t and rj is 
obtained by integrating of Eq. (10) and is given by 23] 


t(i],r) = / a(x, r) dx + h(r) , 


(14) 


where tf,(r) is a functional constant of integration called bang function, since it corresponds to 
the fact that in these models the scalar factor can vanish at different times at different locations. 
We will consider models with tb(r) = 0. The solution of eq.(13) can then be written in the form 

Po 


aM k(r) + 3p(f| 92 (r), 93 (r))’ 
where p(x] g2, g$) is the Weierstrass elliptic function and 

92W = jfc(f) 2 , 9a(r) = (2 k(rf - APo) ■ 


(15) 


(16) 


In terms of r/ and a(r/, r) the radial null geodesics and the Hubble parameter are given by 12] 

dp d T t(r], r) + G(r], r ) 


dz 

dr 


(1 + z)d v G(r],r) 
a{v, r) 


dz (1 + z)d v G(r], r) ’ 
H(rj,r) = H(t(rj,r),r), 


where 


G(r;,r) = 




[5 r (rt(p, r)r) — a 1 d v (a(7},r)r)d r t(i},r)] 


(17) 

(18) 

(19) 

( 20 ) 
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The function G(r), r ) has an explicit analytical form, making the above geodesics equations 
particularly suitable for a low-red-shift expansion. 


4 
















III. LOW RED-SHIFT EXPANSION OF THE HUBBLE PARAMETER H(z) 


In order to obtain a low-redshift formula for the Hubble parameter we expand the function 
k{r ) as 

k(r) = k 0 + kir + k 2 r 2 + ..., (21) 

We also need an expansion for t(rj, r) which can be obtained from the exact solution for a(rj, r ) 
according to 

t{v, r) = t 0 (r) + a(rj 0 , r)(rj - p 0 ) + ^d v a(rj 0 , r){rj - i] a ) 2 + ..., (22) 

where we defined the function t 0 (r ) by 


«o(r) = t(i)o, r). 


(23) 


Using the properties of the Weierstrass elliptic functions p, ( and cr 
integral in eq.(14), obtaining 


42] we can compute the 


t{v,r) = 


2po 
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where p' is the derivative of the Weierstrass’ elliptic function p. p 
the equations: 



( and cr are defined by 


P 1 (p( x )) = x , 


C ( x ) 

o’ (x) 
a(x) 


-fW , 
CW ■ 


(25) 

(26) 

(27) 


The use of the exact expression for t(rj,r) improves the accuracy for the expansion for the 
geodesics respect to previous calculations [24], which were based on a pertubative expansion of 
to(r), rather than the use of the exact value. 

Now we can find the low red-shift Taylor expansion for the geodesic equations 19], and then 
calculate the Hubble parameter. We expand the solution of the geodesic equations according 
to 


r(z) = r\Z + r 2 z 2 + r 3 z 3 + ... 

(28) 

Vi z ) = Vo + Viz + rj 2 z 2 + ... 

(29) 


5 











After substituting the above expansion in the geodesic equations we can map the solution of 
the system of differential equations into the solution of a system of algebraic equations for 
the coefficients of the expansion. Here we give the formulae for the case in which ko = 0, 
which is enough to understand qualitatively the effects of the inhomogeneity. The term k 0 
corresponds in fact to the homogeneous component of the curvature function, which in absence 
of inhomogeneities is simply the curvature of a FLRW model, and as such is not associated to 
any new physical effect not already known from standard cosmology. 

For the geodesics we get: 


Vi 

m 

r\ 

r-2 

r 3 


%H 0 


[°o + *o 


[aoH o t' o (0) (30m (90m — 4) — 8 Ai) + c^Hq (90 2 M — 4Ab) + 


12a%H$Q M 

+ 60m (3 Hq (Om — 1) ^o(9) 2 — ^o(O))], 
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aoHo 
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12<2n-Hf)0 


gUZoliM 
1 


-[a 0 (4AT, - 901,) + 6 (2 - 30 M ) 0 M *o(0)], 


72a 3 0 H 2 Q A n 2 M 


[a 2 0 H 0 (4 K 2 (20a + Co (2 — 30 m) + Om) — 60AiOaO1,+ 


+30a0m ( 8 A 2 + 3 (90m — 4) O^,)) — 36cloHoQ a Qm70'(0) (K\ (40m — 2) + 

+ (8 - 90m) Oi,) + 180 A 0l, (6H q (301, - 40 M + l) t' 0 (0 ) 2 + (2 - 30 M ) %(0 ))], 


where Om, Oa, T 0 and K n are dimensionless quantities given by [16] 


(30) 


(31) 

(32) 

(33) 


(34) 


Po — SQmUqHq , (35) 

A = 30a#o, (36) 

To = Vo (a 0 H 0 ) , (37) 

K n = k n (a 0 H 0 ) n+2 . (38) 


We have also used the following definitions 


a 0 = a(v 0 ,0), 


(39) 


H 0 = H(rj 0 , 0), 

Co = c(y;^.4(2/q-27Sl 4 Jl'|,)) , 


(40) 

(41) 


where ( is the Weierstrass Zeta Function. As we can see in the above formulae the effects of 
the inhomogeneity start to show respectively at first order for rj(z) and second order for r(z). 
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FIG. 1: We plot ^(0) as a function of K\. This is the quantity determining the leading order effect 
for ^ as shown in ea. (l5TT) . 

In order to obtain a formula for the Hubble parameter as a function of the red-shift we need 
to substitute eq. (17-18) in eq. (HUD 

H(z) = H( V (z),r(z)). (42) 

After expanding up to second order in z we get: 

H{z) = H 0 + H\z + H 2 z 2 , (43) 

^ = ^^(^- + 3 ) , (44) 

z V a o J 

H 2 = vo 2 n n—[«o^o (20 (Co — 1) -^1 + 48AiHa^m + 2712a (4 — 3 hi a/) &m) + 

1 za 0 iZAiiM 

+ 6a 0 i7o^A^M^o(0) (20Ai + 9 (8 — 5 11 m) 12 m ) + 1812a12m (7/q (25+ 

-12fl M ) to(0) 2 + 5t"(0))]. (45) 

The procedure to reduce the analytical formula to this form is rather complicated since it 

physically 


42] and we 


involves to express wherever possible all the intermediate expressions in terms o 
meaningful quantities using the properties of the Weierestrass elliptic functions 
give more details in appendix A. 

As we can see from the first order coefficient Hi, at leading order in 2q(0) determines the 
sign of the correction respect to the homogeneous case, and for this reason we plot 2q(0) as 


7 










<5plPb 



FIG. 2: The density contrast 5 = -£■ is plotted as a function of the redshift for three different 
inhomogeneities profiles modeled by LTB solutions. 


a function of K\ in fig.flU). At second order we have a more complicated dependency for H 2 , 
which involves also K 2 and 

We can easily interpret the linear behavior shown fig. 1, applying the chain rule for the 
derivative 


m = 


a = 


dt 0 (r ) dk 


dk 

(a^Ho 


dr r 

-3 ^0 ( r ) 

1 8k 


=0 — otKi , 


\k=ko 


-0.57. 


(46) 

(47) 


IV. RELATIVE DIFFERENCE OF H(z) RESPECT TO THE HOMOGENEOUS 
CASE 


For a flat FLRW solution the expansion rate is given by 

H FLRW {z ) = H 0 y/n M (l + z) 3 + n A . (48) 


Since we want to compare the inhomogeneous case with the flat FLRW case we define the 
relative difference as 


6H(z) _ H ALTB (z) 
~H “ H FLRW (z ) 
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(49) 









where we are denoting with H ALTB (z ) the expansion rate defined in eq. (1421) . 


We can now expand the above expression at low-redshift to get 


SH(z) 5^ 5H 2 2 

H ~ ~ir z + ~W Z + "' 

5H\ 2 i q (0) 2 crOjv/ 


ao 


ao 


K u 


H 

5H 2 _ 

H 36a 2 0 H 0 tt A n M ' 

+9 (8 - 9 n M ) n M ) + 9fi A n 2 M (h 0 (25 - i2U M ) ^(o) 2 + 5t"(o))]. 


(50) 

(51) 


-[2o,qH q Ki (5 (<^o — 1) Af + 1214^12^) + 3czo-^o^a^m7o(0) (20Ai+ 


(52) 


It is straightforward to verify that in the homogeneous limit, when K 0 = K\ = K 2 = 0 , as 
expected, 5Hi = 5H 2 = 0, since £q( 0) = tg(0) = 0 . From the expression for it is clear the 
crucial role played by fg(0), which determines the sign of the relative difference at leading order 
in red shift, and according to eq.(46-47) is proportional to K\. 

Using eq. (S3) we can also see that for a fixed AT 0 the sign of K\ determines the sign of 
at very low red-shift when the second order contributions can be neglected. This in in 
good agreement with fig.([3ll4j), which show an approximate linear behavior with opposite slope, 
corresponding to different values of K\. 


A. Testing the accuracy of the formulae 

In order to test the accuracy of the analytical formulae and compare it with perturbation 
theory we perform numerical calculations using LTB solutions corresponding to the density 
contrasts shown in fig. (El) . Using large density contrasts we can test the limit of the perturbation 
theory results and compare them with the low red-shift expansion. 

As seen in fig.((3jl4j) the analytical formula for the relative difference of H(z) respect to the 
homogeneous case is in good agreement with SH H obtained by integrating numerically the 
geodesics and background equations and is more accurate than the perturbation theory. From 
fig- (SI) in particular we can see that the agreement with the exact numerical calculation for the 
red-shift expansion is in general better than that of the perturbation theory result. For larger 
density contrasts, as expected, the pertubative calculation is increasingly inaccurate while the 
red-shit expansion is still in good agreement with the exact numerical result, since it does not 
rely on the assumption of a small density contrast. This is due to the fact that perturbation 
theory is based on the assumption that 5p/p -C 1, while the low red-shift expansion is based 
on an exact solution of the Einstein’s equations. 
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FIG. 3: The relative difference with respect to the homogeneous case is plotted as a function 

of the redshift for three different density contrasts. The colors correspond to the density contrasts 
in fig.®. The solid lines are for the numerical calculation, the dashed lines are for the analytical 
formulae and the dot-dashed lines are for the perturbation theory result. 

V. CONCLUSIONS 

We have derived a low-redshift analytical formula for the Hubble parameter for an observer 
at the center of a spherically symmetric matter distribution, using an exact solution of Einstein’s 
field equations. Such a formula is in good agreement with numerical calculations and is more 
accurate than the perturbation theory result, especially for large density contrasts. This is due 
to the fact that perturbation theory is based on the assumption that 5p/p <C 1, while the low 
red-shift expansion is based on an exact solution of the Einstein’s equations. 

If the H(z) observations will show deviations from the A CDM predictions compatible with 
the formulae we have derived, this could be considered an independent evidence of the existence 
of a local inhomogeneity, and the formulae could be used to determine the characteristics of 
this local structure. 

While the expansion for r(z ) depends on our choice of radial coordinate, the formulae for the 
H(z) are independent of it, since both H and the redshift are physically observable quantities, 
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FIG. 4: The relative percentual difference A(z) = 100 SH ^ 5H ^ S ^ n ul. of 5H/H respect to the nu¬ 
merical results is plotted as a function of the redshift for the analytical formula (dashed) and for 
perturbation theory (dot-dashed). The colors correspond to the density contrasts in fig.©. As it can 
be seen the low red-shift formula is always better, especially for larger density contrasts. 

and as such are independent of the gauge choice. Since in the derivation of the formulae the 
inhomogeneity profile is determined by the coefficients of the expansion of the function K(r), 
there is still some coordinate dependency left in the way we parameterize the inhomogene¬ 
ity. While our choice of gauge is quite natural, since in the FRW gauge, corresponding to 
Po(r) — const , the radial coordinates reduces to the radial FRW comoving coordinate in the 
limit in which the function k(r) goes to zero, a totally coordinate independent formula would 
still be preferable. For this reason in the future it will be interesting to derive formulae which 
are completely independent of the choice of the coordinate system, parameterizing the inhomo¬ 
geneity in terms of the redshift expansion of the density p(z), which is a physical observable, 
instead of using the expansion of K{r). 
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Appendix A: Derivation of the analytical formula 


In order to obtain the formula for the red-shift expansion of H(z) we have applied several 
manipulations and substitutions. The method is based on re-expressing everything in terms of 
physical quantities, starting from the definitions of a 0 and H 0 , which are related to p and p' 
by the equations 


a o = (Vo,0) = 


H 0 = H(rjo, 0 ) 


P o 

ko + 3po ’ 
- 

2 p 0 ’ 


( 1 ) 

( 2 ) 

(3) 


where 


P o = P('^o;^2(0),^ 3 (0)) , 
dp(v,g2(0),g 3 (0)) 


P o = 


dr) 


\v=vo ■ 


By inverting the previous equations we obtain the following relations 


(4) 

(5) 


Po = p(w,92(0),g 3 (0)) = Po a ° k ° , 

OCLq 

/ dp(r);g 2 (0),g 3 (0)) t 2H 0 p 0 

Po- drj \n=vo~ 3 

We can then substitute the above expressions everywhere p and p' appear, making the final 
formula only depending on physical quantities such as Ho. In order to simplify the results we 
have also used the Einstein’s equation at the center (■ rjo , 0) 


( 6 ) 

(7) 


1 — —Ao + + klj\ . 


( 8 ) 
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